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We discuss the renormalization of the lattice formulation of the Heavy Quark Effective Theory (LHQET). 
In addition to wave function and composite operator renormalizations, on the lattice the classical velocity is also 
renormalized. The origin of this renormalization is the reduction of Lorentz (or 0(4)) invariance to (hyper)cubic 
invariance. We present results of a new, direct lattice simulation of this finite renormalization, and compare the 
results to the perturbative (one loop) result. The simulation results are obtained with the use of a variationally 
optimized heavy-light meson operator, using an ensemble of lattices provided by the Fermilab ACP-MAPS 
collaboration. 



1. INTRODUCTION 

In the lattice formulation of the Heavy Quark 
Effective Theory (LHQET), a new renormalization 
beyond those of the continuum theory arises. It is 
the the renormalization of the "classical velocity". In 
the continuum, the classical velocity v is a fixed 
parameter which appears in the decomposition of 
the momentum of a heavy particle and the reduced 
Dirac equation of the heavy quark field. 

P = Mv + p 

(1) 

-ivDh >v) (x) = 

In the continuum, the second relation is derived from 
the first, and the velocity that appears in these two 
contexts is the same. However, on the lattice this is 
not the case. As we have explained elsewhere, the 
origin of the difference between the values of v in 
these two contexts is the reduction of continuum 
Lorentz (or Euclidean 0(4)) invariance to 
hypercubic invariance on the lattice. 1 

The evaluation of this renormalization is necessary 
in order to relate calculations performed in the 
LHQET to physical quantities. One of the most 
basic of such quantities is the slope of the 
Isgur- Wise form factor, £(vv') at the origin. In the 
M Q -> oo limit, £ describes all semi-leptonic decays 



of heavy mesons into one another, such as the 
process B — > D * / V , The renormalization of v is 
directly reflected in the value of this slope. 

In this talk, we will describe a non-perturbative 
calculation of the leading multiplicative 
renormalization of v by means of a lattice simulation. 
We will first outline how the renormalized value of 
the classical velocity can be extracted from a lattice 
simulation, and then describe a variational procedure 
for enhancing the quality of the signal one obtains in 
the simulation. Finally, we will give the results, and 
compare them to the results of from a one-loop 
perturbative calculation. It is useful to develop the 
calculation in a power series in the "bare" classical 
velocity, and we follow this procedure. 

2. THE PHYSICAL CLASSICAL VELOCITY 

The heavy quark effective theory is efficiently 
formulated 2 by factoring the M — > singular behavior 
from the field of a heavy quark, leaving a reduced 
operator 

h M ( X ) = e 1 y v \|f(jc) (2) 

There is an independent reduced field for each value 
of the classical velocity. To define the classical 



velocity non-perturbatively without global gauge 
fixing, it is useful to consider a gauge invariant 
quantity containing the reduced heavy quark field, 
for example a composite operator for a heavy meson 
made of one light and one heavy quark. The 
asymptotic behavior of its propagator will be 



(3) 



where its rate of fall-off, relative to the 
p = , v = heavy quark energy, is given by 



lim \j (M + m) 2 + ((M + m) v (phys) + p) 2 
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(phys) -(phys) 

= mv + v -p 



Its first derivative with respect to the residual 
momentum is the physical classical velocity 
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3. EXPANSION IN POWERS OF v 

The discretization of the HQET reduced Dirac 
equation we use to put the limiting theory on the 
lattice is 

v [ U(x,x+t)S(x+t,y)-S(x,y) ] 
3 -iv 



[U(x,x+(i) 5(x+p,y) 
z 

-U(x,x-(i) S(x-(i,y) ] 
= §(x,y) 



(6) 



The use of an asymmetric forward time difference 
facilitates implementing the requirement that heavy 
quarks propagate only forward in time. This 
descretization has the virtue that the heavy quark 
propagator is obtained by forward recursion. We 
have described the properties of this implementation 
of the HQET elsewhere. 3 



The bare classical velocity plays the role of a 
transverse hopping constant. There is one factor of v. 

for each step in the i th direction. After after n time 
steps the propagator is an (n - l) st order polynomial 
in the bare classical velocity components. We 
exploit this structure by developing all quantities that 
depend on the classical velocity in a power series in 
the bare components. For example, the meson 
propagator is given as 

A <l) (f,/?) = Y, v i v 2 V HUp,m) (7) 

m 1 ,m 2 ,nt 3 

The computation of the coefficients in this 
polynomial is highly efficient. The indices m. give 

the maximum lattice displacement of the heavy quark 
propagator in the i th direction. 

4. THE OPTIMIZED MESON FIELD 

While in principle any gauge invariant operator 
containing the reduced heavy quark propagator could 
be used to extract the physical calssical velocity, in 
actual simulations it is crucial that the propagator 
reach its asymptotic form in as few time steps as 
possible. We therefore choose a composite field 
which in the Coulomb gauge has the form 4 

*P <v) (x) = X) V{y)q{x)h {v Xx+y) (8) 

y 

and determine the weighting function by the 
requirement that the propagator of this field be 
maximal on the n th time slice. 

The space-time propagator of this composite field 

is 

< ] (jc,jt') =£V(?) (9) 

where 

K^\x,y;x,y') - 



s(x,x)S (v \x +y,x +y) 



(10) 



The average is over the ensemble of lattices. It is 
amusing to note that while we could "center" the 
meson field at the heavy quark, the light quark, or 
anywhere in between, it is handiest to center it at the 



light quark. The requirement that the propagator be 
maximal on a given time slice implies that \\f is that 
eigenvector of with the largest eigenvalue. The 
solution to this eigenvalue problem is greatly 
facilitated by first solving the static, p = , v = 
problem, and then sucessively calculating the 
coefficients of higher powers of v. . 

5. SIMULATION OF v <phys> 

We have computed the magnitude of the classical 
velocity renormalization by simulating Eq. (9) on a 
set of 24 3 x 48 lattices with P = 6.1 and K = .154 
generously provided by the Fermilab Collaboration. 5 
By expanding the logarithmic derivative of the 
optimized heavy meson residual propagator in 
powers of the bare classical velocity, we find the 
combination of terms which should grow linearly 
with lattice time, with slope equal to the negative of 
the coefficient of the first power of the bare classical 
velocity in the physical classical velocity. Taking 
the finite lattice approximation to the momentum 
derivative, the appropriate expression is 



A(t,p.=p,m. = 1)- A(t,p.= -p,m.= 1) 
2p A(t,p=0,m=0) 



(11) 



where p is the minimum momentum on the finite 
lattice. Figure 1 is the graph of this expression, with 
its slope, as fitted either over the full 6 time slice 
range of apparent linearity, or only over the 4 
interior slices. 

It is interesting to compare the result of this 
simulation with the shift in the classical velocity 
evaluated to 1 loop in lattice perturbation theory. 6 
For the forward time-symmetric space discretization 
we have used 



8v. 



.104 + .024 (5 
.121 + .045 (6 

.2335 



10 fit) 
9 fit) 

(one loop) 



„ (bare) 



One-loop lattice perturbation theory grossly 
overestimates the renormalization of the classical 
velocity. 
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Figure 1 — The expression Eq. (11), and its slope 
which gives the leading physical classical velocity 
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